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The dissociative adsorption of hydrogen on Pd(100) has 
been studied by ab initio quantum dynamics and ab initio 
molecular dynamics calculations. Treating all hydrogen de- 
grees of freedom as dynamical coordinates implies a high di- 
mensionality and requires statistical averages over thousands 
of trajectories. An efficient and accurate treatment of such 
extensive statistics is achieved in two steps: In a first step 
we evaluate the ab initio potential energy surface (PES) and 
determine an analytical representation. Then, in an indepen- 
dent second step dynamical calculations are performed on the 
analytical representation of the PES. Thus the dissociation 
dynamics is investigated without any crucial assumption ex- 
cept for the Born-Oppenheimer approximation which is any- 
how employed when density-functional theory calculations are 
performed. The ab initio molecular dynamics is compared 
to detailed quantum dynamical calculations on exactly the 
same ab initio PES. The occurence of quantum oscillations 
in the sticking probability as a function of kinetic energy is 
addressed. They turn out to be very sensitive to the symme- 
try of the initial conditions. At low kinetic energies sticking 
is dominated by the steering effect which is illustrated us- 
ing classical trajectories. The steering effects depends on the 
kinetic energy, but not on the mass of the molecules. Zero- 
point effects lead to strong differences between quantum and 
classical calculations of the sticking probability. The depen- 
dence of the sticking probability on the angle of incidence is 
analysed; it is found to be in good agreement with experi- 
mental data. The results show that the determination of the 
potential energy surface combined with high-dimensional dy- 
namical calculations, in which all relevant degrees of freedon 
are taken into account, leads to a detailed understanding of 
the dissociation dynamics of hydrogen at a transition metal 
surface. 
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I. INTRODUCTION 

The dissociative adsorption of molecules on surfaces is 
one of the fundamental reaction steps occuring in cataly- 
sis, corrosion, and in the hydrogen gas storage in metals. 
Adsorption corresponds to a process in which statistically 
distributed molecules hit the surface from the gas phase. 
For a diatomic molecule this requires the calculation of 
thousands of trajectories; the adsorption probability is 
then obtained by averaging over these trajectories. In 
"traditional" ab initio molecular dynamics the electronic 



structure, total energy and the forces acting on the nuclei 
are determined for each configuration during the journey 
of the particles which requires of the order of 100 - 1000 
self-consistent calculations for each trajectory. Since the 
determination of total energies and forces is still a heavy 
computational task, the number of trajectories obtain- 
able in such a "during the journey" ab initio molecular 
dynamics simulation is limited to numbers well below 
lOO.Ero These numbers are too low to achieve a reason- 
able statistics. Typically one needs to consider of the 
order of 10 3 to 10 5 trajectories, depending on the kind 
of experiment to be simulated.u Therefore we propose a 
"divide and conquer" approach for ab initio molecular 
dynamics: At first the road map should be created and 
only then the journey started. Thus, in a hrst (elab- 
orate) study we analyse the ab initio potential energy 
surface (PES). This PES is parametrized in a suitable 
form and only then the molecular dynamics calculations 
are performed on the analytical representation of the ab 
initio potential. In this way it is easy to study up to 10 6 
trajectories. 

If hydrogen is involved in the dissociation process, 
quantum effects may play a role in the dynamics. We 
have recently improved a quantum dynamical algorithm 
that makes it possible to treat all six degrees of free- 
dom of the hydrogen molecule in the dissociation process 
quantum dynamically.!] Quantum effects can thus be de- 
termined by comparing the results of the quantum dy- 
namical calculations with classical trajectory calculation 
on exactly the same PES. It turns out that for hydrogen 
moving on a strongly corrugated and anisotropic PES 
zero-point effects can be substantial!?] One advantage of 
the quantum dynamics, that to our knowledge has not 
been emphasized in the literature yet, is that the aver- 
aging over the initial conditions is done automatically. 
For example, a plane wave in a j = rotational state 
describing the incident molecular beam hits the surface 
everywhere in the surface unit cell and contains all molec- 
ular orientations with equal probability. 

Dissociative adsorpli^ ..systems can be roughly di- 
vided into two classesElHfl'aEI: Systems, where the stick- 
ing probability is monotonously increasing as a func- 
tion of the incident kinetic energy of the imping- 
ing molecules, and systems, where the sticking prob- 
ability shows an initial decrease with increasing ki- 

first class. These systems are characterized by a min- 
imum barrier hindering dissociative adsorption, so that 
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increasing the kinetic energy helps to overcome the bar- 
rier. 



transition metal surfac- 



dsorption-systems like 
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which is the subject of 
, An initially decreas- 



system H 2 /Pd(100) 
our study, belongs to this class 
ing sticking probability had usually been explained by a 
precursor mechanismn In this concept the molecules are 
temporarily trapped in a molecular physisorption state, 
the so-called precursor state, before they dissociatively 
adsorb. The energy dependence of the sticking probabil- 
ity is related to the trapping probability into the precur- 
sor state. It is this trapping probability which decreases 
with increasing energy.U 

However, it has for example been shown for the system 
H 2 /W(100)-c(2x2)CuE] that for a hydrogen molecule 
impinging on a metal substrate the energy transfer to 
substrate phonons is much too small to account for 
the high sticking probabilities at low kinetic energies 
due to the large mass mismatch. Therefore also di- 
rect non-activated adsorption together with a steering 
effect has been suggested in order to explain the ini- 
tial decrease of the. sticking probability by King almost 
twenty years agce9. Still, in low-dimensional dynami- 
cal treatments of ti^jJt a 7Pd(100) system no steering ef- 
fect was observed Hj'c3'E3 Only very recently it has been 
shown by high-dimensional quantum dynamical calcula- 
tions based on ab initio potential-energy surfaces (EE SI 
for the systems H 2 /Pd(100)l3ffl and H 2 /W(100)BEl 
that indeed an initial decrease of the sticking proba- 
bility with increasing kinetic energy is not necessarily 
due to a precursor mechanism. For both systems the 
PES has non-activated paths towards dissociative ad- 
sorption and no molecular adsorption well. However, 
the majority of pathways towards dissociative adsorption 
has in fact energy barriers with a rather broad distribu- 
tion of heights and positions, i.e. the PES is strongly 
anisotropic and corrugated. Similiar features of the po- 
tential have recently— been found also for the interaction 
of H 2 with Rh(100)B A slow molecule moving on such a 
PES with an unfavorable initial configuration for disso- 
ciative adsorption can be steered efficiently towards non- 
activated paths to adsorption by the forces acting upon 
the molecule. This mechanism becomes less efficient at 
higher kinetic energies because then the molecule is too 
fast to be diverted significantly. More particles are there- 
fore scattered back into the gas-phase from the repulsive 
part of the potential. This leads to the initial decrease 
of the sticking probability. At even higher kinetic ener- 
gies the molecule will eventually have enough energy to 
directly traverse the barriers causing the sticking proba- 
bility to rise again. 

In our calculations for the interaction of H 2 /Pd(100) 
all six degrees of freedom of the hydrogen molecule are 
treated dynamicallya This makes it possible to inves- 
tigate the influence of all hydrogen degrees of freedom 
on the dissociative adsorption, scattering and associative 



desorption on an equal footing. So far we have studied 
the dependence of adsorption and desorptipn on kinetic 
energy, molecular rotation and orienta,tioncrc3, molecular 
vibrationo, ro-vibrational couplingcll and the rotation- 
ally elastic and inelastic diffraction of H 2 /Pd(100)c3. 

In this article we will first describe the quantum and 
classical methods we have used to determine the adsorp- 
tion dynamics of hydrogen on Pd(100). Then we address 
the origin of oscillations in the sticking probability as a 
function of the kinetic energy. Next the steering effect 
is illustrated and the differences between classical and 
quantum dynamics are discussed as are isotope effects in 
the dissociative adsorption dynamics. Finally we will fo- 
cus on the dependence of the sticking probability on the 
angle of incidence. 



II. COMPUTATIONAL DETAILS 

In our approach the dynamical simulations including 
all relevant degrees of freeedom are performed on an 
analytical representation of the ab initio PES. Thus in 
principle we apply only one approximation, namely the 
Born-Oppenheimer approximation, i.e., we assume that 
the electrons follow the motion of the nuclei adiabati- 
cally. Obviously in practice there is a second important 
approximation, namely the treatment of the exchange 
and correlations effects in the density functional theory 
calculations. 

As far as the number of relevant degrees of freedom is 
concerned, in the case of hydrogen dissociation on densely 
packed metal surfaces usually no significant surface rear- 
rangement upon adsorption occurs, and there is only a 
small energy transfer from the light hydrogen molecule 
to the heavy substrate atoms. Even if there is any sur- 
face relaxation upon hydrogen adsorption, it occurs typ- 
ically on a much larger time scale than the adsorption 
event. The crucial process in the dissociative adsorp- 
tion for these particular systems is therefore the conver- 
sion of translational and internal energy of the hydro- 
gen molecule into translational and vibrational energy of 
the adsorbed hydrogen atoms. Thus the dissociation dy- 
namics can be described by a six-dimensional PES which 
takes only the molecular degrees of freedom into account. 
In the following we present our formalism in such a six- 
dimensional formulation. In principle, however, it can be 
extended to include also the substrate degrees of freedom 
if they are relevant. 



A. Parametrization of the ab initio potential 

In order to obtain a reasonable analytical representa- 
tion of the PES, first a sufficient number of ab initio total 
energies has to be computed. High-symmetry points of 
the multi-dimensional configuration space are reflected 
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by extrema in the PES. Typically in between the high- 
symmetry points the PES is smooth and has no addi- 
tional maxima and minima. Of course, this assumption 
has to checked carefully. In the case of a rigid surface, 
the PES of a diatomic molecule interacting with this sur- 
face is a function of the six molecular degrees of freedom. 
They can be represented, e.g., by the center-of-mass co- 
ordinates X, Y, Z , the interatomic distance r, and the 
polar and azimuthal angle of the molecular axis and <j>. 
Figure [l] shows a two-dimensional cut through the six- 
dimensional coordinate space of H 2 /Pd (100), a so-called 
elbow plot. The two considered coordinates are the H 2 
center-of-mass distance from the surface Z and the H-H 
interatomic distance r. 

In order to solve the time-independent Schrodinger 
equation describing dissociative adsorption it is advan- 
tageous to transform the coordinates is. [the-, Zr-plane 
into reaction path coordinates s and pEJ'ES-Eil. Here s 
describes the position along the "reaction path" - the 
dashed line in Fig. |l| - and p is the coordinate perpendic- 
ular to s (see section |IIB| ). We have then parametrized 
the function V(X, Y, s,p, 6, (j>), which describes the po- 
tential energy surface on which the hydrogen molecule 
moves, in the following form 13 

V (X,Y,a,p,6,<f>) = V cori + V lot + V vih (1) 

with 

2 

V co " = V^} n (s) cos mGX cos nGY, (2) 

m,n— 



V TOt = ^ V^\s) - cos 2 B (cos mGX + cosmGY) 

m=0 

2 i 

+ Y V n ] ( S ) 2 Sin2 9 C0S ^ ( C0S UGX ~ C0S nGY ) ( 3 ) 

n=l 

and 

t/vib = ^2 (s) [p Ap{X ^ s)] 2_ (4) 

G = 2n/a is the length of the basis vectors of the square 
surface reciprocal lattice, a is the nearest neighbor dis- 
tant between Pd atoms and lo(s) is ther-wibrational fre- 
quency. We note that Wiesenekker et alEi have recently 
employed an equivalent analytic representation to de- 
scribe the 6D-PES of H 2 /Cu(100), the only difference 
being that they use cartesian coordinates in the Zr-plane 
instead of reaction path coordinates. 

It turns out that the calculation of total energies for 
250 different configurations is sufficient to determine 
the parameters necessary to describe the functions ap- 
pearing in the potentiaL-parametrization. The ab ini- 
tio energies are obtainedEJ using density functional the- 
ory together with the generalized gradient approximation 
(GGA)t3 for the exchange and correlation functinal and 
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FIG. 1. Contour plot of the PES along a two-dimensional 
cut through the six-dimensional coordinate space of 
H2/Pd(100). The inset shows the orientation of the molec- 
ular axis and the lateral H2 center-of-mass coordinates, i.e. 
the coordinates X, Y, 8, and <j>. The coordinates in the figure 
are the H2 center-of-mass distance from the surface Z and 
the H-H interatomic distance r. The dashed line is the opti- 
mum reaction path. Energies are in eV per H2 molecule. The 
contour spacing is 0.1 eV. 

the full-potential linear augmented plane wave method 
(FP-LAPW) (see Ref. || and references therein). The 

ten functions Vm} n (s) and uj(s) are determined such that 
the difference to the ab initio calculations on the average 
is smaller than 25 meV. 

The introduction of the displacement Ap in the poten- 
tial term V vlb (Eq. [|) takes into account that the location 
of the minimum energy path in the Zr-plane depends on 
the cut through the six-dimensional configuration space 
of H 2 /Pd(100). Ap does not influence the barrier dis- 
tribution, however, it changes the curvature of the min- 
imum energy paths in the Zr-planes. For the calculated 
elbow plots Ap reaches values of up to 0.6 A, but only for 
large separations of the hydrogen atoms, i.e., when the 
molecule is already dissociated. Large values of Ap re- 
quire a large number of vibrational eigenfunctions in the 
expansion of the hydrogen wave function in the coupled 
channel scheme (see below), which makes the calculations 
very time-consuming. Since the large values of Ap occur 
only for large separations of the hydrogen atoms, they 
do not influence the calculated sticking probabilities and 
scattering properties significantly, as we have checked by 
test calculations. We have therefore parametrized the 
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displacement properly only for values of |Ap| < 0.15 A. 



B. Quantum dynamics 

The Hamiltonian of a hydrogen molecule interacting 
with a rigid surface can be written as 



H 



2M 



V 2 



R 



— V 

2(i r 



V(R,r) 



(5) 



R = (ri + r 2 )/2 = {X, Y, Z) and r = (r 2 - T X ) are the 
center-of-mass and relative coordinates of the hydrogen 
molecule, respectively, M = 2m and p, = m/2 are the to- 
tal and reduced mass of the hydrogen molecule, respec- 
tively, where m is the mass of a hydrogen atom. Now 
we write the relative part in spherical coordinates; with 
r = |r| this yields 



H = 



t. 

2M ' " 2(i 
V(R,r,9,cj>) , 



X7 2 - — 
V R ~ 



dr 2 



2d_ 

r dr 



l 2 



(6) 



where L is the angular momentum operator. 

We define r e as the equilibrium bond length for a cer- 
tain molecular configuration, i.e., r e is in general a func- 
tion of the five coordinates X, Y, Z, 9 and 4>. The vibra- 
tional amplitude of diatomic molecules is usually small 
compared to the equilibrium bond length, which means 
that 



-| < r 



(7) 



This allows usE3 to neglect the term in the Hamil- 

o r or 

tonian and to approximate the angular momentum term 

T 2 T 2 

by 75- ~ Then we end up with the following Hamil- 
tonian 



H 



V 



R 



2M 
+ V(R,r, 



ti_ 
2(i 



if 



(8) 



As already mentioned in section II A , it is advanta- 



geous to transform the coordinates in the Zr-plane into 
reaction path coordinates s and p in order to solve the 
time-independent Schrodinger equation describing disso- 
ciative adsorption. To perform the transformation to the 
reaction path coordinates, first Z has to be mass-scaled 
to the coordinate z by 



(9) 



The PES in Fig. [l| is already plotted according to mass- 
scaled coordinates. With the reaction path coordinates 
s and p the Hamiltonian becomes 



H 



2(i 
~2p 



V 



d_ 



dp 



'd P/ 



dr 2 

( d 2 



Lr 



2M \dX 2 



d 2 \ 
dY 2 



+ v{x,Y, s , P ,e,4>) ■ (10) 



The coupling parameter rj is defined by 
77 = 1 - k{s) p , 



(11) 



where k(s) is the curvature of the lowest energy reaction 
path (the dashed line in Fig. [I]). 

The displacement Ap also enters the denominator of 
the angular momentum term in the Hamiltonian via the 
equilibrium bondlength r e 



r e = r° e {s) + sin0 r (s) Ap(X,Y,s) 



(12) 



Here sin</v(s) is the angle between the reaction path 
(dashed line in Fig. [I]) and the z-axis. Since the relation 



sin </v(s) Ap{X,Y,s) < r° e {s) 



(13) 



2 in the angular momentum 
(s) Ap(X, Y, s) up to sec- 



holds, we have expanded r t 
term in a Taylor series in sin 
ond order. 

The quantum dynamical calculations are performed by 
solving the time-independent Schrodinger equation for 
the two hydrogen nuclei moving on the six-dimensional 
PES in a coupled-channel scheme. As channels the eigen- 
functions of the Hamiltonian for molecules in the gas 
phase are used. We.-use.the concept of the local reflec- 
tion matrix (LORE)E3C3. For a detailed description of 
this stable coupled channel method we refer to Refs. ^ 
and In the LORE scheme the reflection matrix R is 
determined; in order to obtain sticking probabilities Si 
for some initial state i, where i stands for a multi- index, 
we use unitarity: 



St = 1 - 



(14) 



Rji is the differential reflection amplitude; the sum over 
j extends over all possible reflection states. 

The basis set used in the coupled-channel algorithm for 
the H2 results presented here included rotational eigen- 
functions with rotational quantum numbers up to j max = 
8, vibrational eigenfunctions with vibrational quantum 
numbers up to Umax = 2, and parallel momentum states 
with maximum parallel momentum p max = 7hG with 
G = 2-Kja. The convergence of the results with respect 
to the basis set has been carefully checked by calculations 
with maximum quantum numbers j max — 10, v max = 3, 
and p max = 10?iG, respectively. 

Due to the higher mass of D2 the energy spacing be- 
tween the quantum levels is smaller for D2 than for 
H2. Therefore much more eigenfunctions in the expan- 
sion of the wavefunction have to be taken into account 
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in the coupled-channel calculations for D2 than for H2. 
This makes a six-dimensional quantum treatment of D2 
very time-consuming. In order to investigate isotope ef- 
fects we have therefore performed five-dimensional vibra- 
tionally adiabatic quantum calculations for D2 , where the 
molecules are not allowed to make vibrational transitions. 
We have already shown that vibrationally adiabatic cal- 
culations are very close to the full six^dimensional results 
for the dissociation of H2 on Pd(100)EJ. This should also 
be valid for D2 since the ratio of the vibrational time- 
scale to the other time-scales of rotation and translation 
is the same for H2 and D2 . The five-dimensional quantum 
calculations for D2 have been performed with rotational 
quantum states up to j max = 12 and parallel momentum 
states with p max = WTiG. 

In coupled-channel calculations always the whole S- 
matrix has to be computed. This leads to a ./V 3 -scaling 
of the algorithm due to the matrix operations, where N 
is the number of channels included in the calculation. 
In order to perform these demanding quantum dynam- 
ical calculations it is therefore necessary to utilize the 



symmetries of the scattering problem (see also Ref. 57). 



First of all selections rules are important. Because of the 
inversion symmetry of the H2 molecule only rotational 
transition with Aj — even are allowed, where j is the 
rotational quantum number. In addition, the analytic 
representation of PES only contains rotational potential 
terms that cause Am = even transitions (see the second 
sum of V rot in Eq. ^), where m is the azimuthal quantum 
number of the H 2 molecule. 

Furthermore, we exploit the symmetry group of the 
Hamiltonian which corresponds to the Ci v symmetry of 
the fee (lOO)-surface. In general the scattering solutions 
do not belong to irreducible representations of the sym- 
metry group. However, if the scattering solutions of inter- 
est can be decomposed into irreducible representations, 
the number of relevant channels per coupled-channel cal- 
culation can be significantly reduced. This is due to 
the fact that only channels belonging to the same ir- 
reducible representation of the symmetry group couple 
to each other, since the Hamiltonian commutes with the 
symmetry operator. 

If, for example, the incident parallel momentum corre- 
sponds to a reciprocal lattice vector (this includes the 
zero-vector for normal incidence) and the initial rota- 
tional quantum number j and the azimuthal quantum 
number m are even, the scattering solutions can be bro- 
ken up into eight different irreducible representations of 
the symmetry group, four of which can be further de- 
composed. In each decomposition the number of chan- 
nels is roughly halved, and in each irreducible representa- 
tion the S'-matrix is calculated separately. This leads to 
a reduction of the computational cost to approximately 
4 • (1/8) 3 + 8 • (1/16) 3 = 5/512 w 1.0%. If only the 
sticking probability for normal incidence is required, it is 
sufficient to calculate only two S- matrices, i.e., the ex- 
ploitation of the symmetries causes a reduction in the 
CPU time to - 2 • (1/16) 3 = 1/2048 « 0.05%. With- 



out the use of the symmetry the calculations presented 
here would not be feasible. Using the selection rules and 
the decomposition into irreducible representations up to 
25,000 channels per total energy are taken into account in 
the quantum dynamical calculations; the actual number 
of channels in the single calculations is usually < 600. 
For a more detailed description of the construction of 
symmetry adapted channels, see Ref. E7L 



C. Classical dynamics 

The classical trajectory calculations are performed on 
exactly the same PES as the quantum dynamical calcu- 
lations. To derive the classical equations of motion isoni 
the reaction path Hamiltonian Eq. we have useefcia 



-ih d s = p s , 
-ifi dp=p p . 



(15) 



The equations of motions are numerically integrated with 
the Bulirsch-Stoer method with a variable time-stepEa. 
We required that the energy conservation per molecu- 
lar dynamics run was fulfilled to 0.1 meV. The sticking 
probability is determined by averaging over a sufficient 
number of trajectories. The exact number of trajectories 
to be considered depends on the specific initial conditions 
and ranges between 1,815 and 18,330. 

As far as the CPU time requirement is concerned, it 
is a wide-spread believe that classical methods are much 
less time-consuming than quantum ones. This is cer- 
tainly true if one compares the computational cost of one 
trajectory to a quantum calculation. However, since in 
quantum mechanics the averaging over initial conditions 
is done automatically while in classical mechanics one has 
to average over many trajectories corresponding to dif- 
ferent initial conditions, for the results presented here 
the quantum method is even more time-efficient than 
the classical calculations, in particular if one considers 
the fact, that in a coupled-channel method the sticking 
and scattering probabilites of all open channels are de- 
termined simultaneously. 



III. RESULTS AND DISCUSSION 
A. Quantum oscillations 

Figure |^ presents six-dimensional quantum dynamical 
calculations of the sticking probability as a function of 
the kinetic energy of a H2 beam under normal incidence 
on a Pd(100) surface and five-dimensional calculations 
for D2. In addition, the results of the H2 molecular beam 
experiment by Rendulic, Anger and WinkleiO are shown. 

First of all a very strong oscillatory structure is ap- 
parent in the sticking probability as a function of the 
incident energy. Such structures reflect the quantum na- 
ture of the scattering. They are known for a long time in 
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FIG. 2. Sticking probability versus kinetic energy for a 
hydrogen beam under normal incidence on a Pd(100) sur- 
face. Experiment (H2): circles (from ref. theory: 
six-dimensional results for H2 molecules initially in the rota- 
tional and vibrational ground state (dashed line) and with an 
initial rotational and energy distribution adequate for molecu- 
lar beam experiments (solid line), and vibrationally adiabatic 
five-dimensional results for D2 molecules initially in the rota- 
tional ground state (dash-dotted line). 

He and H2 scattering^ and also in LEED0. In the case 
of the sticking probability of H2/Pd(100), these -i^adlla.- 
tions have been the issue, of a current debateHCd'EjO. 
We have recently shownO that the peaks in the sticking 
probability can be related to the opening up of new scat- 
tering channels with increasing kinetic energy, especially 
at low kinetic energies. In particular the emergence of 
the [10], [11], and [20] diffraction channels for normal in- 
cidence and the opening up of rotational excitation lead 
to strong peaks in the sticking probablity. Here [n, m] 
corresponds to the diffraciinu. indices of the (lOO)-surface. 

Rettner and AuerbachEdo have tried to find the the- 
oretically predicted oscillationsQ by an effusive beam ex- 
periment, but they could not detect any. As pointed 
out^al£3 the height of the peaks is very sensitive to the 
symmetry of the scattering conditions. Surface imperfec- 
tions like adatoms and steps and also the thermal motion 
of the substrate will reduce the coherence of the scatter- 
ing process and thus smooth out the oscillatory structure. 
But more importantly, also the angle of incidence has a 
decisive influence on the symmetry 

The experiment by Rettner and Auerbach was done 
for an angle of incidence of 15°, while the calculations 
were done for normal incidence. In order to investigate 
the dependence of the oscillatory structure on the angle 
of incidence we performed five-dimensional vibrationally 
adiabatic calculations for an angle of incidence of the 
molecular beam of 15° and compared them with normal- 
incidence results (see Fig. ph. The energy resolution for 
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FIG. 3. 5D quantum calculations of the sticking probability 
versus kinetic energy for a hydrogen beam impinging on a 
Pd(100) surface. Upper panel: H2 under normal incindence 
0i = 0°; Middle panel: H2 with an incident angle of 9i = 15°; 
Lower panel: D2 under normal incidence 0i = 0° . 



the non-normal incidence results was chosen to be below 
the width of the prominent peaks for normal incidence 
so that these peaks should be detected. In order to rule 
out basis set effects, also the results for normal incidence 
in Fig. U were obtained by five-dimensional calculations. 
Caused by the reduced dimensionality, the height of the 
peaks for normal incidence is changed compared to the 
full six-dimensional calculations. This indicates the im- 
portance of the full dimensionality for the calculations. 
The peak positions, however, are the same. Also the av- 
eraged sticking probability is not changed significantly^. 

Now at normal incidence every diffraction channel is at 
least fourfold degenerate (except for the backscattered 
beam) due to the C± v symmetry of the (100) surface. 
This makes the effect of the opening up of a new scat- 
tering channel much more dramatic than in the case of 
a general angle of incidence, where this degeneracy is 
lost. This is demonstrated in Fig. ||, which shows that 
in the energy regime belo w 40 . meV, which was probed 
by Rettner and Auerbach,L3o the sticking probability 
for 9i = 15° is much smoother than for normal inci- 
dence. Thus, as expected, calculations for normal inci- 
dence are not directly comparable with experiments for 
non-normal incidence, at least as far as quantum effects 
are concerned. 

The large peak at approximately E.- L w 50 meV is 
still visible for both angles of incidence. This peak is 
due to the opening up of rotationally inelastic diffrac- 
tion, i.e., the kinetic energy becomes large enough to en- 
able j = — > 2 rotational transitions in scattering. For 
9{ = 15° this peak is slightly shifted to higher total ki- 
netic energies. This is simply due to the fact that the 
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rotationally inelastic specular peak opens up at higher 
total kinetic energy due to the parallel momentum con- 
servation. 

In recent quantum dynamical calculations of the dis-. 
sociative adsorption of the reactive system H2/W(100)e3 
also oscillations in the sticking probability were found, 
but they were much smaller than for H2/Pd(100). This 
may be caused by the lower dimensionality of these cal- 
culations. Since only one surface coordinate was con- 
sidered, the number of degenerate scattering channels 
opening up was less than in six-dimensional calculations, 
leading to smaller effects. For example, in previous quan- 
tum dynamical calculations of the dissociative adsorp- 
tion of H2 on a model potential with activated as well as 
non-activated paths to adsorption, where also only one 
surface coordinate was considered, quantum oscillations 
have also been foundjj there amplitude, however, is much 
smaller than in the 6D-calculations. 

For the heavier isotope D 2 the energetic spacing be- 
tween quantum levels is much less due to the higher mass 
compared to H2 . The higher "density" of channels should 
also make the effects of the opening up of new scattering 
channels less dramatic. In the lower panel of Fig. || we 
have plotted five-dimensionals results for the D2 sticking 
probability under normal incidence with the same energy 
resolution as for the H2 non- normal incidence results. Ex- 
cept for low kinetics energies the sticking curve is indeed 
much smoother than for H2 at normal incidence. 

There is a further source for the smoothening of the 
sticking probability as a function of the incident kinetic 
energy in supersonic molecular beam experiments: The 
experimental beam does not correspond to a monoener- 
getic beam in one specific quantum state. Instead, there 
is a certain velocity spread of the impinging molecules 
which is typically of the order of Av/vi =0.1, where Vi 
is the mean initial velocityeS; in addition, the internal 
states of the molecules are populated according to some 
Boltzmann-like distribution. For the solid line in Fig. ^ 
we have assumed an initial rotational and energy distri- 
bution adequate for molecular beam experiments. As a 
consequence, the oscillatory structure is almost entirely 
washed out. Accordingly, also the experimental data of 
Ref. ^5| do not show any significant oscillations. 

B. Steering effect 

We will now discuss the general trends in the aver- 
aged sticking probability as a function of the kinetic en- 
ergy. The qualitative features of the experimental stick- 
ing probabilityc3 are well reproduced by the averaged 
quantum dynamical results, as Fig. || shows, although 
there are quantitative differences which we will address 
below. At low energies there is a substantial decrease 
in the sticking probability with increasing kinetic energy, 
which is then followed by an increase at higher kinetic 
energies. As already pointed out above, such a gen- 




FIG. 4. Snapshots of classical trajectories of hydrogen 
molecules impinging on a Pd(100) surface. The initial con- 
ditions are chosen in such a way that the trajectories are re- 
stricted to the zz-plane. Left trajectory: initial kinetic energy 
Ei = O.Of eV. Right trajectory: same initial conditions as in 
the left trajectory except that the molecule has a higher ki- 
netic energy of 0.12 eV. 

eral behavior had usually been attributed to a precursor 
mechanism, in which the impinging molecules prior to 
dissociation are first trapped in a physisorption well due 
to energy transfer to substrate phononsH. This mecha- 
nism, however, cannot explain the quantum dynamical 
results since first, there is no physisorption well in the 
calculated PES, and second, there is no energy transfer 
to the surface possible due to the use of a fixed substrate. 
Thus the decrease in the sticking probability has to be 
caused.iby.-a uurplv dynamical effect, namely the steering 
effect :Bc3u'OCjcj Although the majority of pathways to 
dissociative adsorption has non- vanishing barriers with a 
rather broad distribution of heights and positions, slow 
molecules can be very efficiently steered to non-activated 
pathways towards dissociative adsorption by the attrac- 
tive forces of the potential. This mechanism becomes less 
effective at higher kinetic energies where the molecules 
are too fast to be focused into favourable configurations 
towards dissociative adsorption. If the kinetic energy 
is further increased, the molecules will eventually have 
enough energy to directly traverse the barrier region lead- 
ing to the final rise in the sticking probability. 

In order to illustrate the steering effects, we use the re- 
sults of two typical classical trajectory runs. This is done 
in Fig. |^, where snapshots of these two trajectories are 
shown. The initial conditions are chosen in such a way 
that the trajectories are restricted to the xz-plane. The 
left trajectory illustrates the steering effect. The incident 
kinetic energy is Ei = 0.01 eV. Initially the molecular 
axis is almost perpendicular to the surface. In such a 
configuration the molecule cannot dissociate at the sur- 
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face. But the molecule is so slow that the forces acting 
upon it can reorient the molecule. It is turned paral- 
lel to the surface and then follows a non-activated path 
towards dissociative adsorption. 

In the case of the right trajectory, the initial condi- 
tions are the same as in the left one, except that the 
molecule has a higher kinetic energy of 0.12 eV. Due to 
the anisotropy of the PES the molecule also starts to ro- 
tate to a configuration parallel to the surface. However, 
now the molecule is so fast that it hits the repulsive wall 
of the potential before it is in a favorable configuration 
to dissociative adsorption. It is then scattered back into 
the gas-phase rotationally excited. 

Fig. H shows that there are still quantitative differences 
betweeen theory and experiment. Considering the fact 
that there are no adjustable parameters in our calcula- 
tions, the agreement is quite satisfactory, though. The 
discrepancies might be due to uncertainties in the deter- 
mination-of the ab initio PES which are of the order of 
±0.1 eVEil. We also like to point out that the experi- 
mental values rfithe sticking probability are subject of a 
current debateEj'Eil. 

Furthermore, in our calculations substrate phonons or 
electronic excitations are not taken into account. We 
have noted above that due to the large mass mismatch 
between impinging hydrogen molecule and the Pd sub- 
strate atoms the substrate motion can be neglected as 
far as understanding the basic dissociation mechanism 
is concerned. Taking the substrate motion into account 
would allow for recoil of the surface atoms upon impact of 
the impinging molecules. Although the energy transfer 
to the solid is rather small, recoil of the surface atoms 
leads to_a_sticking curve which is stretched to higher 
energiesEJO. In other words, it renormalizes the energy 
axis, because due to the energy transfer to the surface the 
effective kinetic energy becomes smaller. Such a renor- 
malization would improve the agreement between exper- 
iment and theory, as an inspection of Fig. || reveals. 



C. Comparison quantum-classical dynamics and 
isotope effects 

In Fig. [5] we compare the averaged quantum mechani- 
cal sticking probability for H2 with the results of classical 
and quasiclassical trajectory calculations for H2 and D2. 
The inset shows an enlargement of the results at low ener- 
gies. Quasiclassical in this context corresponds to trajec- 
tories with the initial vibrational energy of the hydrogen 
molecule equal to the vibrational zero-point energy of hy- 
drogen, which is 0.258 eV for H 2 and 0.185 eV for D 2 , 
while for the purely classical trajectories the molecules 
are initially non-vibrating. First of all, the classical re- 
sults do not show any oscillatory structure revealing that 
the oscillations are entirely due to quantum mechanics. 
Note that the quasiclassical calculations for H2 show al- 
most no initial decrease in the sticking probability. For 




c 

o 

H — • 

CO 



0.2 



0.0 



0.00 0.05 0.10 

— 6D quantum dynamics H 2 , beam simulation 
O Quasiclassical H 2 
♦ Classical H 2 
■ Quasiclassical D 2 
A Classical D, 



0.0 0.2 0.4 0.6 0.8 

Kinetic energy E ; (eV) 



1.0 



FIG. 5. Probability for dissociative adsorption versus ki- 
netic translational energy for a hydrogen beam under nor- 
mal incidence on a clean Pd(100) surface for non-rotating 
molecules. The solid line shows six-dimensional quantum 
dynamical results for H2 assuming an energy spread typical 
for beam experiments. The quasiclassical results corresponds 
to molecules initially vibrating with an energy equal to the 
vibrational zero-point energy while the classical results are 
obtained with initially non-vibrating molecules. Open dia- 
monds: quasiclassical H2; filled diamonds: classical H2; filled 
squares: quasiclassical D2; open triangles: classical D2. The 
inset shows an enlargement of the results at low energies. 



D2 there is a small decrease, while the purely classical 
results effectively fall upon the averaged H2 quantum re- 
sults at low and high kinetic energies. 

We have recently shown that the strong difference be- 
tween quasiclassical results on the one side and classi- 
cal and averaged quantum results on the other side is 
caused by zero-point effects of the hydrogen-. molecule 
in the multi-dimensional configuration spacetl. When 
the molecule approaches the surface, the molecular bond 
is weakened and consequently the molecular vibration 
is softened, i.e., the vibrational frequency decreases. 
Since the change of the frequency is slow compared 
to the vibrational period, the vibrational energy pfcd-. 
lows the change of the frequency almost adiabaticallyL]'E3 
which leads to an effective vibrational-translational en- 
ergy transfer. At the same time, due to the anisotropy 
and corrugation of the PES the molecule about to dis- 
sociate becomes localized in the remaining four degrees 
of freedom of the molecule which are the polar and az- 
imuthal rotation and the two translations parallel to the 
surface. This localization leads to the building up of 
additional zero-point energies due to the Heisenberg un- 
certainty principle. In fact, the sum of all zero-point en- 
ergies remains approximately constant alons-a-4ninimum 
energy path towards dissociative adsorptiont3E3, and for 
H 2 /Pd(100) the sum becomes even larger than the gas- 
phase vibrational zero-point energy of H2, which is the 
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only zero-point energy of a free molecule. 

Now a classical particle can follow precisely a mini- 
mum energy path through a corrugated PES; in a picto- 
rial sense one might say it can propagate along the bot- 
tom of the valley in the PES. A quantum particle cannot 
do that. It has to be delocalized and needs to have at 
least the zero-point energies perpendicular to the prop- 
agation direction to traverse a corrugated PES without 
tunneling. This leads to an effective upwards shift of the 
potential for the quantum particle along the minimum 
energy path. Or vice versa, the classical particle experi- 
ences a lower minimum potential. As a consequence, the 
sticking probability for the quasiclassical particle is much 
larger than for the quantum particle. 

As mentioned above, the sum of all zero-point energies 
along the minimum energy path towards dissociation ad- 
sorption increases in the systen H2/Pd(100), but in a first 
approximation we assume it to be roughly constant. In 
such a situation the combined effect of all zero-point en- 
ergies is to cause a constant shift of the potential. There- 
fore the results of quantum calculations and classical cal- 
culations without any initial zero-point energy should be 
similiar since a constant shift of the potential does not 
affect the dynamical properties. This is indeed the case 
at low and at high energies, as the comparison between 
purely classical and quantum results in Fig. demon- 
strates. In addition, these results confirm that steering 
is a general dynamical effect and is not particularly re- 
lated to quantum or classical dynamics. 

The problem of a proper treatment of zero-point en- 
ergies in quasiclassical trajectory calculations, is. well- 
known, especially in the gas-phase communitjoo One 
possible way to deal with this problem is the reduced 
dimensionality treatment in the vibrationally adiabatic 
approximation (for a overview see Ref. 71). In this ap- 
proach a small number of degrees of freedom is treated 
dynamically while the remaining degrees of freedom are 
taken into account by adding the sum of their zero-point 
energies to the potential along the reaction path. An- 
other more elaborate approach is to constrain the energy 
in each, idbrational mode to be greater than its zero-point 

In our purely classical approach we ignore zero-point 
energies all along the reaction path. But this approach 
is actually in the spirit of the vibrationally adiabatic ap- 
proximation. It effectively takes the zero-point energies 
into account through a shift of the potential along the 
reaction path corresponding to the sum of all zero-point 
energies. This shift, however, is constant along the reac- 
tion path. Moreover, we still keep the full dimensionality 
of the problem by explicitly treating all degrees of free- 
dom dynamically. This is indeed essential since for ex- 
ample the steering effect is absent inL^_Law-dimensional 
treatment of the H 2 /Pd(100) systemBEiH 

Besides zero-point effects tunneling is also an impor- 
tant quantum phenomenon. However, in a system with 
activated and non-activated paths towards dissociative 
adsorption tunneling does not play an important role. 



This is due to the fact that tunneling is exponentially 
suppressed. Hence the propagation of the quantum par- 
ticel along a classically possible path is much more prob- 
able than the dissociation via tunneling. 

The results also show that in purely classical dynam- 
ics there is no isotope effect between H2 and D2 in the 
sticking probability. As far as the low-energy regime is 
concerned, this seems surprising at a first glance, since 
D2 is more inert than H2 due to its higher mass. How- 
ever, one has to keep in mind that at the same kinetic 
energy D2 is slower than H2, so that there is more time 
for the steering forces to redirect thc-.D2 molecule. This 
has been noted before by Kay et allEi. 

Indeed, the Lagrangian for a system of classical parti- 
cles with the same mass Mi can be written as 



L 



EM\ ( dxi 
~Y \~~dT 



U({xi}) 



(16) 



For another isotope with the mass M 2 the potential does 
not change. If we transform the time axis via 



(17) 



we end up with the following Lagrangian for the new 
isotope of mass M2: 



V 



E 



Mi ( dx t 

He 



U({xi}) 



(18) 



which is equivalent to the Lagrangian of Eq. [l^. This 
means that the equations of motion for an system of clas- 
sical particles with mass Mi correspond to the equations 
of motion for a system of classical particles with mass M2, 

where the velocities have been scaled by an factor J 

i.e., where the kinetic energy is the same. Hence, for 
different isotopes with the same initial conditions, where 
only the initial velocities have been scaled to keep the ki- 
netic energy unchanged, the trajectories remain exactly 
the same. 

It follows that there cannot be any isotope effects as a 
function of the kinetic energy for hydrogen moving clas- 
sically on a PES as long as there are no energy transfer 
processes to, e.g., substrate phonons. Furthermore, this 
indicates that the steering effect is not restricted to light 
molecules as hydrogen, but should also be operative for 
all other heavier molecules moving in a similiar PES. As 
far as dissociative adsorption is concerned, however, for 
heavier molecules recoil effects of the substrate might no 
longer be negligible, so that for a complete description 
considering all relevant degrees of freedom energy trans- 
fer processes to substrate phonons could be important. 

Contrary to the purely classical calculations, the qua- 
siclassical results show an isotope effect between H2 and 
D 2 . The sticking probability of H2 is larger compared 
to D 2 , the effect being most pronounced for kinetic ener- 
gies between 0.03 eV and 0.30 eV. This isotope effect can 







only be caused by the different initial vibrational zero- 
point energies which can be effectively used to traverse 
the corrugated and anisotropic barrier region. The H2 
gas-phase zero-point energy is larger by 73 meV; indeed 
the H2 sticking curve seems to be shifted to lower energies 
with regard to the D2 sticking curve by approximately 
this amount. 

Interestingly enough, at very low kinetic energies below 
0.03 eV also the quasiclassical calculations show almost 
no isotope effect, in addition to the fact that classical 
and quasiclassical results are almost identical in this low- 
energy range. In the limit of zero initial kinetic energy 
apparently the sticking probability is to a large extent 
determined by steering forces which already act rather 
far away from the surface where the change of the vi- 
brational frequency and thus zero-point effects are in- 
significant. The absence of an isotope effect for very low 
energies is actually also true for the averaged quantum 
results, as Figs. || and || demonstrate. However, there is a 
pronounced isotope effect in the quantum results for ki- 
netic energies larger than 0.1 eV. The size of this isotope 
effect corresponds to the one found in the quasiclassical 
calculations which again shows that the different initial 
vibrational zero-point energies cause the isotope effect. 

We also like to comment on the rather large differ- 
ence between classical and quantum results in Fig. ^| for 
kinetic energies between 0.15 eV and 0.6 eV. We think 
that this difference might be due to the fact that the 
sum of all zero-point energies along the minimum en- 
ergy path through the barrier region actually becomes 
larger than the gas-phase zero-point energy.cl This effect 
is most prominent in the medium energy range, where 
steering is no longer effective. At very high kinetic ener- 
gies, where zero-point effects should play only a negligi- 
ble role, indeed quantum and classical results are in very 
close agreement. Furthermore, in quantum mechanics it 
always takes a finite amount of energy to change the state 
of a particle (if there are no degenerate states), while in 
classical mechanics particles can be diverted by any in- 
finitesimally small amount of energy. This should make 
the quantum propagation somehow stiffer than classical 
propagation. This could also contribute to the difference 
between quantum and classical results in the medium en- 
ergy range. 



D. Non- normal incidence 

In this section we address the issue of non-normal in- 
cidence. The experiments on the angular dependence 
of the sticking probability of H2/Pd(100) were done for 
two different-jnitial kinetic energies, Ei — 0.1 eV and 
Ei = 0.4 eVEj The incident azimuth was not identified. 
We have determined the quantum and classical angular 
dependence of the sticking probability for two different 
incident azimuths: along the [10] direction which corre- 
sponds to one axis of the surface square lattice, and along 
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FIG. 6. Probability for dissociative adsorption of H2 on 
Pd(100) versus angle of incidence for an inital total kinetic 
energy of Ei =0.1 eV. Experiment: diamonds (from ref. 
The theoretical results are for initially non-rotating molecules. 
Circles show 6D quantum results, squares classical trajectory 
calculations. The filled and open symbols correspond to cal- 
culations with the azimuthal angle of incidence along the [10] 
and [11] direction of the (100) surface, respectively. 



the [11] direction which corresponds to the diagonal of 
the surface square lattice. 

Fig. ^| compares the experiment with the quantum and 
classical results at the lower kinetic energy, Ei = 0.1 eV. 
Note that the quantum results were determined for a mo- 
noenergetic beam in one specific quantum state which is 
here the vibrational and rotational ground state; hen 
quantum oscillations are superimposed on these dat 
but apparently their size is small. The general experi- 
mental trend is well reproduced by the theoretical results: 
the sticking probability is almost independent of the an- 
gle of incidence at this energy, but slightly increases with 
increasing angle. There is no large difference between 
quantum and classical results. Only at angles larger than 
45° the classical results are above the quantum results. 
There is also almost no significant dependence on the 
azimuth except for the classical results at almost graz- 
ing incidence of 9{ = 80° where the sticking probability 
along the [11] direction is larger by 0.1 compared to the 
[10] results. 

The angular dependence of the sticking probability at 
the higher kinetic energy of Ei = 0.4 eV is plotted in 
Fig. 0. First of all the absolut values determined by ex- 
periment, quantum and classical calculations are very dif- 
ferent at this energy, as is already apparent from Fig. ^ 
and Fig. [|; this issue was discussed above. But the 
general trends in the angular dependence are in good 
agreement. All sets of data show a significant decrease 
in the sticking probability with increasing incident an- 
gle for angles below 60°. For almost grazing incidence 
there is now a substantial dependence on the initial az- 
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FIG. 7. Probability for dissociative adsorption of H2 on 
Pd(100) versus angle of incidence for an inital total kinetic 
energy of Ei = 0.4 eV. The meaning of the symbols is the 
same as in Fig. ^[ 

imuth. The sticking probability for molecules impinging 
along the [11] direction is signifincantly larger than for 
molecules impinging along the [10] direction. This differ- 
ence is more pronounced for the quantum than for the 
classical calculations. 

Our results show that the general features of the an- 
gular dependence of the sticking probability determined 
in the experiment - an increase with increasing incident 
angle at low energies and a decrease at higher kinetic 
energies - is well reproduced by our six-dimensional cal- 
culations. In particular, the calculations demonstrate 
that a sticking probability increasing with increasing in- 
cident angle-in not necessarily indicative of a precursor 
mechanismc3 but can be caused by the dynamics of the 
dissociative adsorption on a corrugated PES. 

Still the questions remains: what causes the different 
angular dependence at these two energies? It is useful 
to discuss angular effects by eppsidering the energy scal- 
ing of the sticking probability!?!, i.e., by determining the 
exponent n such that 

S(E i ,6 i )nS(E i cos n 6 i ,6 i = 0°) (19) 

If n — 2, then the so-called normal energy scaling is 
valid, i.e., the sticking probability is a function of the 
normal component of the incident energy alone. In our 
calculations the sticking probability for normal incidence 
has its minimum at approximately Ei « 0.1 eV (see 
Fig. ^) . If normal energy scaling were fulfilled in the sys- 
tem H 2 /Pd(100), then for E, < 0.1 eV the sticking prob- 
ability would indeed rise with increasing incident angle 
since the normal energy decreases, and it would fall with 
increasing angle for Ei > 0.1 eV as long as the normal 
component Ei cos 2 9i is larger than 0.1 eV. 

In order to check whether normal energy scaling is valid 



FIG. 8. Quantum mechanical sticking probability of H2 on 
Pd(100) as a function of the normal component of the incident 
energy for molecules initially in the vibrational and rotational 
ground state. Solid line: beam under normal incidence with 
an energy spread typical for molecular beam experiments; cir- 
cles, diamonds, and squares: monoenergetic beam with an an- 
gle of incidence between 0° and 30° , between 30° and 60° , and 
between 60° and 90° , respectively. The results for non- normal 
incidence are obtained for different azimuthal angles. 

in the system H2/Pd(100), we have plotted the quan- 
tum mechanical sticking probability as a function of the 
normal component of the incident energy Ei cos 2 &i in 
Fig. ||. The non-normal incidence data show some scat- 
ter, in particular at low energies. This can be caused by 
quantum oscillations; the azimuthal dependence, which 
is not specified in Fig. ||, also contributes to the scat- 
tering of the data. But the general trends are in qual- 
itative agreement with model calculations on a three- 
dimensional PES with activated as well as non-activated 
path towards dissociative adsorptiontl At low normal ki- 
netic energies below 0.05 eV additional parallel momen- 
tum suppresses the sticking, for normal kinetic energies 
between 0.05 and 0.35 eV additional parallel momentum 
enhances sticking, and above 0.35 eV the results show 
approximate normal energy scaling. These results also 
show that it is not possible to assign any global energy 
scaling, i.e. any global exponent n in Eq. [l^, to the an- 
gular dependence of the sticking probability. 

At very low energies, where the steering mechanism 
is operative, additional parallel momentum hinders dis- 
sociation. This can be easily understood. A molecule 
with a low normal velocity may still be steered to a fa- 
vorable site for dissociation. But due to the additional 
parallel momentum the molecule is swept past this fa- 
vorable site and scattered back into the gas-phase from a 
repulsive site before the bond-breaking can occur. This 
effect is similiafc 1 to || thc rotational hindering in the steer- 
ing regimeQCj'OOllJ which is caused by the fact that 
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rapidly rotating molecules rotate out of favorable orien- 
tation for dissociation during the interaction with the 
surface. However, Fig. [7] shows that for incident an- 
gles above 70° the suppression of the steering depends 
strongly on the incident azimuth. 

In the intermediate energy range between 0.05 and 
0.35 eV additional parallel momentum enhances sticking, 
in particular for incident angles above 60° . For molecules 
impinging on the surface under an angle larger than 60° 
the component of the kinetic energy parallel to the sur- 
face is at least three times larger than the normal com- 
ponent. These molecules experience an lateral average of 
the PES in this energy range.u Steering in the angular 
degrees of freedom can still occur. Indeed the sticking 
probability for 9{ > 60° shows a decrease for normal ki- 
netic energies between 0.05 eV and 0.12 eV indicating 
a steering effect, and then an increase at higher ener- 
gies. Far away from the surface the molecules are first 
attracted to the on-top-siteC-i But molecules steered to 
this site will eventually encounter a barrier towards disso- 
ciative adsorption of 0.15 eV. In order to dissociate slow 
molecules have to be re-directed towards the bridge or 
hollow sites (see also Ref. [ll]). Thus potential gradients 
can also steer molecules to "wrong" sites. This oversteer- 
ing in the lateral coordinates cannot occur for molecules 
experiencing a laterally averaged potential causing the 
increase in the sticking probability for large additional 
parallel momentum. 

In the direct dissociation regime for normal kinetic en- 
ergies larger than the lateral average of the barrier height 
additional parallel momentum causes an increase in the 
sticking probability!!! This lateral average still depends 
on the orientation of the molecule. For the majority 
of molecular orientations the laterally averaged barrier 
heights for H 2 /Pd(100) are less than 0.15 eV (see, e.g., 
the barrier distribution in Ref. ^). Hence the lateral 
averaging also leads to an increase in the sticking proba- 
bility in the direct dissociation regime. 

This mechanism, however, does not promote stick- 
ing significantly any more if the normal kinetic en- 
ergy is larger than most of the maximum barriers for 
fixed molecular orientation. Still, the fact that for 
Ei cos 2 9i > 0.35 eV the sticking probability shows ap- 
proximate normal energy scaling in spite of the strong 
corrugation of the PES is reminiscent of the activated 
system H2/CU. Thesa-similiar results have been found 
both experimentallytlHi3 and theoretically! 2 !] although the 
PES is also strongly corrugatecOO. This apparent con- 
tradiction is attribJU±ed-.to features of the PES called bal- 
anced corrugationQOEa. For this type of corrugation the 
higher barriers have to be farther away from the surface 
compared to the lower barriers. These features are also 
present in the system H 2 /Pd(lQ0) where the highest bar- 
riers are over the on-top-sitescil. 

We now return to the discussion of the angular de- 
pendence of the sticking probability for fixed total ki- 
netic energy. For fixed total kinetic energy increasing the 
incident angle means decreasing the normal kinetic en- 



ergy and increasing the incident parallel momentum. At 
low kinetic energies decreasing the normal kinetic energy 
makes the steering more effective which promotes dissoci- 
ation. On the other hand, increasing the incident parallel 
momentum hinders dissociation in the low-energy range. 
At Ei = 0.1 eV both effects approximately cancel which 
leads to a sticking probability almost independent of the 
incident angle (see Fig. ||). 

At normal energies larger than 0.1 eV decreasing the 
normal kinetic energy leads to a decrease in the sticking 
probability, but increasing the incident angle enhances 
the sticking probability for normal energies below 0.4 eV. 
However, the promoting effect of additional parallel mo- 
mentum is less pronounced than the decrease due to the 
smaller normal kinetic energy. Hence in Fig. [7] the stick- 
ing probability decreases for increasing incident angle at 
an initial total kinetic energy of 0.4 eV. 

Figs. ^ and ^ also show that for 8i < 60° there is al- 
most no dependence of the sticking probability for non- 
normal incidence on the azimuth. For larger incident 
angles, however, molecules impinging along the [11] di- 
rection of the surface unit cell have a higher dissociation 
probability than molecules impinging along the [10] di- 
rection. This can be explained by a shadowing effect. For 
molecules approaching the surface under an almost graz- 
ing incidence along one axis of the quadratic surface unit 
cell, the most favorable adsorption path at the bridge po- 
sition is effectively hidden behind the high barriers at the 
on-top position. For an approach along the diagonal of 
the square unit cell this most favorable adsorption path 
is still directly accessible. 

Finally we study the influence of the incident rotational 
quantum state on the angular dependence of the sticking 
probability. We have determined the sticking probability 
as a function of the angle of incidence for molecules ini- 
tially in the rotational quantum state jj = 2, m, = and 
ji = 2,771; = 2 for Ei = 0.1 eV. The results are plotted 
in Fig. H Note that for 6>; = 0° the sticking probability 
for initially non-rotating molecules with ji — (Li s (j ■ 3 i (f gc. 
Fig. |). Fig. I shows the well-known resultHMaSEl 
that rotational motion hinders the dissociation at low 
energies because rotating molecules rotate out of fa- 
vorable orientations for dissociation. This suppression, 
however, depends on the orientation of the molecules. 
Molecules with azimuthal quantum number m = j have 
their axis preferentially oriented parallel to the surface. 
These molecules rotating in the so-called helicopter fash- 
ion dissociate more easily than molecules rotating in the 
cartwheel fashion (m = 0) where the rotational axis is 
preferentially oriented parallel to the surface. The latter 
have a high probability hitting the surface in an upright 
orientation in which they cannot dissociate. 

This steric effect is effective for all incident angles, i.e., 
the sticking probability for m t = 2 is always larger than 
for rrii = 0. Like for the non-rotating molecules at this ki- 
netic energy, the results show only a weak dependence on 
the incident angle. For incident angles below 45° there is 
also no significant dependence on the azimuth, but again, 
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FIG. 9. Quantum dynamical probability for dissociative 
adsorption of H2 on Pd(100) versus angle of incidence for 
an inital total kinetic energy of Ei = 0.1 eV. The results are 
for molecules initially in the rotational quantum state j = 2. 
Circles show results for initial rotational azimuthal quantum 
number vm = 0, squares for mi = 2. The filled and open 
symbols correspond to calculations with the azimuthal an- 
gle of incidence along the [10] and [11] direction of the (100) 
surface, respectively. 



for almost grazing incidence molecules approaching along 
the diagonal of the surface unit cell have a higher dissoci- 
ation probability than molecules approaching along one 
axis of the surface unit cell. These results indicate that to 
first order rotational and parallel motion are decoupled 
as far as the dissociation dynamics is concerned. 



IV. CONCLUSIONS 

In conclusion, we reported a six-dimensional quantum 
and classical dynamical study of dissociative adsorption 
of hydrogen on Pd(100). We used a potential energy 
surface obtained from detailed density functional the- 
ory calculations for the system H2/Pd(100). The six 
hydrogen degrees of freedom are treated fully dynam- 
ically. The two main approximations are, firstly, that 
the substrate is kept fixed so that no thermal disorder 
or phonon excitations are allowed, and secondly, that 
the system is assumed to remain in its electronic ground 
state. Hence the continuous excitation spectrum of the 
semi-infinite substrate is neglected. Still these calcula- 
tions are able to reproduce all of the known experimen- 
tal results with regard to the disscociative adsorption at 
least semi-quantitatively. The time-reverse process to 
dissociative adsorption, the associative desorption, was 
not discussed in this study, but previous studies showed 
that also experimental desorption properties are well- 
reproduced by our calculations. Among the processes 



that are now quite well understood are the dependence 
of adsorption and/or desorption on the molecular trans- 
lational energy, vibrational and rotational state, orienta- 
tion of the molecule, and the angle of incidence. 

Quantum effects are non-ncgliblc for hydrogen disso- 
ciation on surfaces. The discrete nature of diffraction 
and molecular excitation leads to a strong oscillatory 
structure of the sticking probability as a function of the 
incident kinetic energy. Furthermore, zero-point effects 
cause substantial deviations between averaged quantum 
dynamical calculations and quasiclassical calculations, in 
which the initial conditions correspond to a molecule vi- 
brating with the gas-phase zero-point energy of hydrogen. 
The corrugated and anisotropic potential energy surface 
leads to the building up of additional zero-point ener- 
gies which effectively increase the minimum potential in 
the quantum calculations. This changes the dynamics in 
the low-energy regime of H 2 /Pd(100) dramatically. How- 
ever, the building up of the additional zero-point energies 
roughly cancels and even overcompensates the decrease 
in the zero-point energy of the H-H vibrations upon dis- 
sociative adsorption. Therefore purely classical calcula- 
tions which neglect the zero-point energies in the initial 
conditions are closer to the quantum results than the 
quasiclassical calculations. 

At low kinetic energies the dissociative adsorption is 
dominated by the steering effect. For higher kinetic en- 
ergies steering becomes less efficient leading to the initial 
decrease in the sticking probability. The steering effect 
is dependent on the kinetic energy, but not on the mass 
of the molecule. Hence steering should also be effective 
for heavier molecules. 

There are still some quantitative differences between 
theory and experiment. They might be caused by uncer- 
tainties in the evaluation of the PES, but also by uncer- 
tainties in the experimental determination of the sticking 
probability. In addition, the differences might be caused 
by the neglect of substrate phonons or electronic excita- 
tions in the calculations. Hence we will address the role 
of the substrate degrees of freedom in the adsorption and 
desorption processes in the future. As for now, our results 
show that the ah initio determination of the potential en- 
ergy surface combined with high-dimensional dynamical 
calculations, in which the relevant degrees of freedon are 
taken into account, is an important step forward in our 
understanding of simple reactions at surfaces. 
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